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We determined the probability distribution of the combined output power from twenty five coupled 
fiber lasers and show that it agrees well with the Tracy-Widom, Majumdar-Vergassola and Vivo- 
Majumdar-Bohigas distributions of the largest eigenvalue of Wishart random matrices with no 
fitting parameters. This was achieved with 500, 000 measurements of the combined output power 
from the fiber lasers, that continuously changes with variations of the fiber lasers lengths. We 
show experimentally that for small deviations of the combined output power over its mean value 
the Tracy-Widom distribution is correct, while for large deviations the Majumdar-Vergassola and 
Vivo-Majumdar-Bohigas distributions are correct. 



Random matrix theory has been exploited in numerous 
research fields ranging from nuclear spectra to quantum 
transport, models of quantum gravity in two dimensions, 
mesoscopic non-linear dynamics, atomic physics, wire- 
less communications and multi-dimensional data analy- 
sis [TH5] . Of special interest are the minimal and maximal 
eigenvalues of random matrices, that determines for ex- 
ample, the conductance fluctuations in two- and three- 
dimensional Anderson insulators [6l [7]- An analytical 
expression describing typical deviations of the maximal 
eigenvalue was presented in the 90's by Tracy and Widom 
(TW) jSI [9] initiating many further theoretical develop- 
ments in random matrix theory [lOl [TT]. Recently, Ma- 
jumdar and Vergassola (MV) calculated the probabil- 
ity of large deviations of the maximal eigenvalue [12- 
[14j above the mean and Vivo, Majumdar and Bohigas 
(VMB) calculated the probability below the mean. The 
MV and the VMB distributions were numerically con- 
firmed, but so far eluded experimental demonstration. 

In this letter, we provide the first experimental obser- 
vation of the MV and VMB distributions in a physical 
system and connect the field of coupled random lasers 
to random matrix theory. We report our measured dis- 
tribution of the combined output power from an array 
of 25 coupled fiber lasers whose cavity lengths randomly 
fluctuate in time. We found that the measured distri- 
bution of the combined output power agrees well with 
the distribution of maximal eigenvalue of Wishart ran- 
dom matrices as predicted by TW and MV. For devia- 
tions close to the mean value, the measured distribution 
is shown to have a universal shape that agrees with the 
TW distribution. For large deviations from the mean 
value the measured distribution deviates from the TW 
distribution, but agrees well with the MV and VMB dis- 
tributions over more than five decades with no fitting 
parameters. To account for this agreement, we present 
a heuristic model that illustrates the relation between 
the output power distribution from our array of coupled 
lasers to the maximal eigenvalue of Wishart random ma- 
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trices. 

Our experiment consisted of an array of 25 coupled 
fiber lasers schematically presented in Fig. [T] Each fiber 
laser was comprised of a Ytterbium doped double clad 
fiber with lengths that varied from 1.3m to 1.7m, a high 
reflecting fiber Bragg grating (FBG) at the rear end of 
the fiber and a low reflecting FBG at the front of the 
fiber. Each fiber lasers was end pumped by a stabilized 
diode laser of 975nm wavelength. The length of each 
fiber laser was about 5m and the output wavelength was 
lOTOnm with a bandwidth of lOnm. Accordingly, there 
were 100,000 available frequencies (longitudinal modes) 
for each laser. The light emerging from all the fiber lasers 
was collected with a lens that was focused onto a detector 
to obtain the combined output power. The fiber lasers 
were arranged in 5 x 5 array, where the coupling between 
them was achieved by means of four coupling mirrors. 
By controlling the orientations of the coupling mirrors we 
could realize a variety of connectivities for the fiber lasers 
in the array, and in our experiments we concentrated on 
the one-dimensional and two-dimensional connectivities. 
Details about the experimental configuration, coupling 
arrangement and connectivity manipulations were pre- 
sented in previous work |15 . 

The lasers were operated close to threshold to maxi- 
mize mode competition and ensure that lasing will only 
occur at the mode where the losses are minimal [161 Hi] • 
We measured the combined output power from the ar- 
ray over a duration of 60 hours. The correlation time of 
the output power fluctuations was found to be shorter 
than 0.5s, hance we obtained over 500, 000 uncorrelated 
measurements. Representative results of the combined 
output power over its mean as a function of time, with 
and without coupling between the lasers, are presented 
m Fig. (2] These are shown over a relatively short time 
duration, but their behavior was similar throughout the 
60 hours measurement. As seen, the power fluctuations 
with coupling (blue curve) are much larger than those 
without coupling between the lasers (red curve), indicat- 
ing that the fluctuations result from the coupling between 
the lasers \1S\ . 

Next we compared the measured results to the dis- 
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where t is the largest eigenvalue and N is the matrix 
size. As evident, there is a very good agreement between 
the probability distributions of the experimental results 
and the TW distribution both for the one-dimensional 
connectivity and for the two-dimensional connectivity. 



FIG. 1. Experimental arrangement for measuring the com- 
bined output power distribution of 25 coupled fiber lasers. 
FBG - fiber Bragg gratings that serve as rear mirror (> 99% 
reflectivity) and front mirrors (- 5% reflectivity); E^^^ - the 
complex electric field in the i'th fiber near the rear FBG for 
each fiber laser; k - the length of the i'th fiber; Mij - cor- 
responds to the propagation matrix for a single round trip 
in the cavity and includes the propagation in each fiber, the 
output coupler 2% reflectivity) and the coupling between 
the different fibers 8% coupling). The light emerging from 
all the fiber lasers was collected with a lens that was focused 
onto a detector to obtain the combined output power. Details 
about the experimental configuration, coupling arrangement 
and connectivity manipulations were presented in previous 
work [15]. 
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FIG. 2. Representative experimental results of the combined 
output power from the 25 fiber lasers as a function of time. 
Red curve - without coupling between the lasers; blue curve - 
with coupling between the lasers. These results indicate that 
the fluctuations are caused by the coupling between the lasers. 



tribution of the largest eigenvalues of Wishart random 
matrices. Figure [3] presents the probability distributions 
of the measured output power in a one-dimensional con- 
nectivity (circles) and in a two-dimensional connectivity 
(asterisks) where the position of the maximum is chosen 
according to the maximum of the TW distribution (solid 
curve) [19 . We present the TW distribution using the 



FIG. 3. Probability distribution of the scaled combined 
output power. Circles - experimental results in a one- 
dimensional connectivity; asterisks - experimental results in 
a two-dimensional connectivity; solid curve - Tracy- Widom 
(TW) distribution. As seen, there is a very good agreement 
between the probability distributions of the experimentally 
measured results and the TW distribution in linear scale for 
both connectivities. For closer inspection of the fit to the tails 
of the distribution see Fig. ^ 

For closer inspection of the tails of the measured dis- 
tributions, we present in Fig. [4] the probability distri- 
butions of the measured combined output power, for 
one-dimensional connectivity (circles, Fig^a)) and two- 
dimensional connectivity (asterisks. Fig. [4p3)), together 
with the TW distribution (solid curves) and MV and 
VMB distributions (dashed curves) on a logarithmic 
scale. The insets illustrate the connectivities between 
the 25 fiber lasers [15]. 

The VMB distribution, plotted with no fitting param- 
eters, is 

P{x) = ci exp [-7V^<I>_ {-x)] when x < 0; (2) 

and the MV is [12] 

P{x) = C2 exp [-N <!>+ {x)] when x > 0; (3) 

where ci = 0.5 and C2 = 0.0063 were found using nu- 
merical simulation in [12 , and the functions ^-\-{x) and 
<l>_(x) are 

M-) = | + l-ln(. + 4) + ^G(^); (4) 
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and 

with G{z) = 3F2[{1, 1, 3/2}, {2, 3}, a hypergeometric 
function. 

As evident from Fig.|4j there are significant systematic 
deviations of the measured distribution from the TW dis- 
tribution, both at values which are much larger or much 
smaller than the mean value p!2HT4] . However, there 
is a very good agreement between the experimental re- 
sults and the MV and the VMB distributions for both 
connectivities, without any fitting parameters. The ex- 
perimental results of the one-dimensional and the two- 
dimensional connectivities are essentially identical indi- 
cating the universality of the maximal eigenvalue distri- 
bution. 

In order to illustrate the relation between the distri- 
bution of the measured power fluctuations and the max- 
imal eigenvalues of Wishart random matrices, we devel- 
oped a simple linear model. While an array of coupled 
fiber lasers is essentially a nonlinear system, many of its 
properties can be determined by its linearized round trip 
propagation matrix [16]. For example, the eigenvectors 
of this matrix correspond to the various global modes of 
the array while the eigenvalues are = 1 — a^, where 
an is the loss of mode n [21 . The tendency of lasers 
to minimize losses will lead the coupled lasers to oper- 
ate in the eigenmode corresponding to the largest eigen- 
value [HI nil. 

We start by letting the electric field near the rear 
FBG of each fiber laser be a component of a vector of 
the total input field |£^o) (see Fig. [T]). After propagat- 
ing one round trip, the field can be described as 
= M l^o), where M is a 25 x 25 round trip propa- 
gation matrix. Details on the derivation of a round trip 
propagation matrix are presented in [21]. Specifically, 
the elements along the diagonal of M denote the self- 
feedback light for each laser, as 

Mi^i = {1 - 4f^)e^'^^^\ (6) 

where hi is the coupling strength between two adjacent 
lasers, li the length of the i'th fiber laser and ki the wave 
vector of the i'th laser out of all the 100,000 available 
frequencies. The off-diagonal elements of M denote the 
coupling between the lasers. For adjacent lasers which 
are not coupled the corresponding elements are zero. 
However, when the adjacent lasers are coupled the cor- 
responding elements above the diagonal are 

Mi^j = i^e'^^^^^^^^\ (7) 

and below the diagonal the elements are 

M^,^- = f^e'^^^^^^^^\ (8) 

In a resonant cavity at steady state, the vector \Ei) 
should be one of the eigenvectors of M, so \Ei) = \Eo) 
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FIG. 4. Probability distribution of the scaled combined out- 
put power for (a) One-dimensional connectivity (circles), and 
(b) two-dimensional connectivity (asterisks) in logarithmic 
scale. Solid curves - Tracy- Widom (TW) distribution; dashed 
curves - Majumdar-Vergassola (MV) distribution for eigenval- 
ues much larger (right green) and Vivo, Majumdar and Bo- 
higas distribution for eigenvalues (VMB) much smaller (left 
blue) than the mean, with no fitting parameters. As seen, at 
the tails of the measured distribution there are significant sys- 
tematic deviations from the TW distribution. However, there 
is a very good agreement between the measured results and 
the MV distribution both at values which are much larger 
or much smaller from the mean value, for both connectivi- 
ties. Insets illustrate the connectivities between the twenty 
five fiber lasers in each case [15]. 



where A^ is inversely proportional to the losses in a sin- 
gle round trip. For many round trips, a complex A^ will 
increase the losses due to interference. These considera- 
tions imply that to ensure minimal losses in the combined 
cavity A^ should be real and maximal [16 . Due to the 
mode competition between the eigenvectors of M on the 
nonlinear gain, the coupled lasers will lase at the mode 
with the minimal losses [17 , which corresponds to the 
eigenvector of the round trip propagation matrix with 
the maximal real eigenvalue. For high gain lasers such as 
fiber lasers, the output power of a mode is proportional 
to its eigenvalue [21]. Therefore, the combined output 
power of the array provides a direct measure for the value 
of the largest eigenvalue of the propagation matrix. 
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Due to thermal and acoustical fluctuations the length 
of each flber laser changes rapidly such that kki mod 
27r is effectively a random phase [181 ES]. These ran- 
dom phases, after each variation in the flber lengths, re- 
sult in a different random round trip propagation matrix. 
The time scale of the length fluctuations in our system 
is much longer than the relaxation oscillation time of the 
lasers [16 , justifying the steady sate assumption. Ac- 
cordingly, the distribution of the combined output power 
from the array flts the distribution of the largest eigen- 
value of random matrices. 

The probability for flnding a single common wave vec- 
tor k such that all the lasers in the array will have the 
same phase and a real valued is negligible — expo- 
nentially proportional to the number of lasers and is 10~^ 
for 25 lasers [24l[25]. Instead, the lasers group in several 
clusters where each cluster has its own wave vector [26] . 
Since the light that is coupled from one cluster to the 
other is essentially lost, the structure of the round trip 
propagation matrix M is block diagonal, where the ele- 
ments along the diagonal are 1 — 4^ and the off diago- 
nal elements when there is coupling between two speciflc 
lasers are So after each fluctuation we have a differ- 
ent matrix where the blocks sizes and locations and the 
signs of the off diagonal elements are random. To show 
that such a round trip propagation matrix M fall on the 
Wishart random matrix universality class, we simulated 
10^ different random realizations of our array with small 
(~ 10/im) fluctuations in the lengths of each flber. In 
each realization we found the clusters with common wave 
vector that yield minimal losses and obtained the corre- 
sponding round trip propagation matrix M [2T. Since 
M represents the round trip propagation in the cavity we 
can deflne a matrix X which represents a single pass in 
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the cavity, so, M = XX^. We evaluated the probability 
distribution of each element in X and found it to have 
a Gaussian shape (data not shown), indicating that the 
round trip propagation matrix M is indeed a Wishart 
matrix. Therefore, the distribution of the combined out- 
put power from the array should flt to the TW and MV 
distributions. 

To conclude, we presented an experimental conflgu- 
ration of 25 coupled flber lasers and showed that the 
probability distribution of their combined output power 
agrees well with the distribution of the largest eigenvalue 
of Wishart random matrices, namely the Tracy- Widom, 
Majumdar- Vergassola and Vivo-Majumdar-Bohigas dis- 
tributions. We believe that such a conflguration can be 
extended to investigate symplectic and non-Hermitian 
random matrices with various connectivities by varying 
the polarizations and the losses in the flber lasers. More- 
over, while in this letter we investigated the combined 
output power from an array of coupled lasers operating 
close to threshold, it is possible to operate the lasers far 
above their threshold and to investigate the phase lock- 
ing across the array [151 El] • Such measurement of phase 
locking gives a direct measure for the number of lasers in 
each cluster and thereby enables investigation of coupled 
ensembles of oscillators where a common frequency for 
all the oscillators can not be found. 
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